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Abstract 

In a variety of supersymmetric extensions of the Standard Model, the scalar part- 
ners of the quarks and leptons are predicted to be very heavy and beyond the reach of 
next-generation colliders. For instance, the realization of electroweak baryogenesis in 
supersymmetry requires new sources of CP-violation, which can only be naturally ac- 
commodated with electric dipole moment constraints if the first and second generation 
scalar fermions are beyond the TeV scale. Also in focus-point supersymmetry and split 
supersymmetry the scalar fermions are very heavy. In this work, the phenomenology 
of scenarios with electroweak baryogenesis and in the focus point region at the LHC 
and ILC is studied, which becomes challenging due to the presence of heavy scalar 
fermions. Implications for the analysis of baryogenesis and dark matter are deduced. 
It is found that precision measurements of superpartner properties allow an accurate 
determination of the dark matter relic density in both scenarios, while important but 
only incomplete information about the baryogenesis mechanism can be obtained. 



1 Introduction 



Super symmetry is one of the most compelling extensions of the Standard Model, with the 
possibility to explain the stabilization of the electroweak symmetry breaking scale, the exis- 
tence of dark matter in the universe and the generation of the baryon-antibaryon-asymmetry 
at the electroweak phase transition (electroweak baryogenesis) . The existence of dark matter 
in the universe has been firmly established by various experiments, and has been measured 
precisely by the Wilkinson Microwave Anisotropy Probe (WMAP) [1], in agreement with the 
Sloan Digital Sky Survey (SDSS) [2], ficDM^ 2 = 0.1106^;^ at the 68% C.L. Here O CD m 
is the ratio of the dark matter energy density to the critical density p c = 3Hq / (87tGn) , 
where H Q = h x 100 km/s/Mpc is the Hubble constant and is Newton's constant. In 
the Minimal Supersymmetric extension of the Standard Model (MSSM), the stable lightest 
neutralino is an attractive candidate for dark matter. 

The process of electroweak baryogenesis requires a sufficiently strongly first order elec- 
troweak phase transition, v(T c )/T c > 1, where v (T c ) denotes the Higgs vacuum expectation 

value at the critical temperature T c [3]. Moreover, new sources of CP violation in addition 
to the CKM matrix phase are necessary [4]. Loop effects of light scalar top quarks (stops) 
can induce a strongly first order electroweak phase transition, thus generating the out-of- 
equilibrium condition for electroweak baryogenesis [5-8]. In addition, CP violation in the 
chargino sector of the MSSM can explain the magnitude of the baryon asymmetry. 

The parameter space for successful electroweak baryogenesis in the MSSM is already 
highly constrained. The lightest stop mass needs to be in the range of roughly 115 to 140 
GeV, whereas the Higgs boson involved in the electroweak phase transition must be lighter 
than about 120 GeV [5-9]. Furthermore, a CP- violating phase in the chargino sector is 
highly constrained by bounds on electric dipole moments [10, 11]. 

Most MSSM scenarios predict a dark matter density larger than the measured value. 
However, the requirement of a light stop from the baryogenesis mechanism suggests that co- 
annihilation between the stop and the lightest neutralino can bring down the relic density to 
the proper region. It has been shown that the co-annihilation is efficient for mass differences 
between the light stop and the lightest neutralino, that are smaller than about 30 GeV [12]. 

Light stops can be searched for at the Tevatron, with a reach up to stop masses of about 
170 GeV for 2-4 fb" 1 of integrated luminosity [13]. However, if the stop-neutralino mass 
difference is less than 30-50 GeV, the stop signal cannot be identified at the Tevatron, due 
to the reliance on a trigger for missing transverse energy. At the next generation of colliders, 
the situation looks more promising. The Large Hadron Collider (LHC), depending on the 
region of parameter space, can see a signal from stops in gluino cascades, while a future 
international e + e~ linear collider (ILC) has excellent capabilities to discover and analyze 
light stops [9,14-16]. 

In order to allow sufficiently large CP violating phases for baryogenesis in the chargino 
sector, but evade current experimental electric dipole moment bounds [17], without invoking 
miraculous cancellations, the sfermions of the first two generations need to be heavier than 
a few TeV, thus effectively decoupling from collider experiments. This situation has grave 
consequences for the investigation of supersymmetry at the LHC and ILC. While the lightest 
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Higgs boson, and possibly some of the heavier Higgs boson, are within experimental reach as 
in other MSSM scenarios, only few supersymmetric particles are likely to be kinematically 
accessible in this scenario: neutralinos and charginos, the light stop and potentially the 
gluino and some sbottom and stau states. Thus it becomes much more difficult to identify 
and measure the properties of these particles. A similar situation occurs in focus-point [18] 
and split supersymmetry scenarios [19], where all sfermions are predicted to be very heavy. 

In this work, it is investigated how a scenario with heavy sleptons and squarks can be 
studied at the LHC and the ILC The analysis is carried out in detail for an electroweak 
baryogenesis scenario and for a focus-point scenario. Split supersymmetry is not studied 
separately, but the conclusions are very similar to the focus-point case. After reviewing the 
definition of the relevant parameters as well as experimental and theoretical constraints in 
section [21 the phenomenology of the baryogenesis scenario is discussed in section [3J while 
section El is devoted to the focus-point scenario. The results are based on a phenomenolog- 
ical analysis, including backgrounds and systematic experimental limitations, but without 
a realistic detector simulation and based on tree-level formulae only. For both scenarios, 
the cosmological implications to be gained from collider results are derived for elucidating 
the nature of baryogenesis and dark matter. For the main part of the work, it is assumed 
that both electron and positron beam polarization are available at the ILC. In section [5] it 
is studied how the results change without positron polarization. Finally, the conclusions are 
given in sectional 

2 Baryogenesis and focus-point scenarios in the MSSM 

Baryogenesis in the MSSM requires an additional source of CP violation beyond the 
Standard Model CKM matrix. Within the MSSM, the dominant source are chargino and 
neutralino loops, with a contribution proportional to Im{/iM 12 } [20,21]. Here \x is the 
Higgs/higgsino parameter and M 2 and M\ are the soft SU(2) and U(l) gaugino parameters, 
respectively. To generate a sufficiently large baryon asymmetry, the charginos are required 
to be relatively light, m-± ~ C(a few 100 GeV). In addition, the CP- violating phase needs 
to be sizable, arg(/iM 2 ) > 0.1 [20]. 

A very large CP- violating phase, on the other hand, is restricted by experimental bounds 
on the electric dipole moments of the electron, neutron and 199 Hg nucleus [17]. The leading 
contributions from one-loop sfermion-gaugino loops [10, 11] become small for large masses 
of the first two generation sfermions of several TeV or for large cancellations between the 
sfermion mass and mixing parameters. Here, to avoid the constraints from the one-loop 
contributions, the sleptons and squarks of the first two generations are assumed to have 
masses of about 10 TeV. Then the leading contributions from supersymmetric particles to the 
electric dipole moments come from two-loop diagrams. For moderate values of tan (3 < 10, 
the ratio of the vacuum expectation values of the two Higgs doublets, and the CP-violating 
phase, it has been shown [22, 23] that the two-loop effects are in agreement with current 
experimental constraints. 

Generically, at the LHC many supersymmetric particles could be observed in the decay 
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chain of squarks, whereas at ILC selectron production would be one of the cleanest and 
most precise testing grounds for supersymmetry. However, in the baryogenesis motivated 
scenario, the large mass of the first and second generation sfermions effectively decouples 
them from observable processes. This has a large impact on the prospects for future collider 
experiments, since the squark and slepton production channels mentioned above are not 
available. Likewise, the dynamics of baryogenesis and dark matter annihilation are governed 
by the light particles in the MSSM spectrum, which include, beside the stops, the charginos 
and neutralinos. Therefore the following analysis focuses on the phenomenology of the 
chargino and neutralino particles. 

The spectrum of the two charginos and four neutralinos is described by the parameters 
p, Mi, M 2 and tan/?. A CP- violating phase in M 2 can always be transferred into the p 
parameter by means of a unitary transformation. In principle, there can also be non-trivial 
phases in the gaugino parameters Mi and M 3 . While the effect of a phase of M 3 on elec- 
troweak baryogenesis is small, a complex phase of Mi could have interesting consequences, 
but is not investigated further here. Therefore in the following all gaugino soft parameters 
are assumed real, while the generation of the baryon asymmetry is connected with a phase 
in the p parameter, p = \p\ x e"^\ 

For definiteness, the specific MSSM parameter point BGEN will be considered in the 
following, as defined in the appendix. Since the allowed parameter space for baryogenesis in 
the MSSM is already highly constrained by experimental bounds, this particular scenario is 
typical for the general MSSM baryogenesis case. 

In the MSSM Higgs sector, the tree-level masses of the CP-even neutral Higgs bosons 
h° and H° and the charged scalar H^ 1 can be expressed through the gauge boson masses, 
the mass of the pseudo-scalar Higgs boson, m^o , and tan (3. The Born relations are however 
significantly modified by radiative corrections, with dominant effects originating from top 
and stop loops. The mass of the Higgs state related to electroweak symmetry breaking (in 
most cases h°) is very sensitive to the stop spectrum. In order to be consistent with the 
bound m-ifi > 114.4 GeV from direct searches at LEP [24] and with one light stop state, 
the heavier stop mass has to be above about 1 TeV and the trilinear coupling A t has to be 
sizable [7]. Constraints from electroweak precision data, in particular the p parameter, are 
satisfied when the light stop is mainly right-chiral. This is naturally achieved for values of 
the stop supersymmetry breaking parameters > 1 TeV 2 and < 0, respectively. The 
stop mixing parameter X t = p cot (3 — A t is bounded from below by the Higgs boson mass 
constraint from LEP and from above by the requirement of the strength of the first order 
electroweak phase transition, leaving the allowed range 0.3 < \X t \/m^ < 0.5 [7]. The value 
of tan/3 is also constrained to the range 5 < tan/3 < 10, with the upper bound stemming 
from present electric dipole moment limits and the lower bound is again related to the LEP 
Higgs mass limit. The latter can be weakened for large values of wIq 3 of several TeV. Also 
values of m^o larger than about 200 GeV are preferred in order to be compatible with the 
electric dipole moment bounds. 

The MSSM Higgs masses with CP violation have been calculated including complete one- 
loop and leading two- loop corrections, see e.g. Ref. [25]. In this work, however, the process of 
baryogenesis at the electroweak phase transition is computed with the program of Refs. [20, 
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22], which includes only one-loop corrections to the zero temperature Higgs potential. Since 
the allowed mass range for the Higgs boson is constrained by the mechanism of electroweak 
baryogenesis, for consistency the Higgs mass is determined by the minimization of the one- 
loop effective potential. This implies that only one-loop corrections are included in the 
calculation of the Higgs mass as well. Although when including two-loop corrections, the 
values of the fundamental parameters will change for given values of m^o and v(T c )/T c , the 
correlation between m h o and v(T c )/T c is expected not to be strongly modified 1 . 

Focus point supersymmetry [18] was suggested to solve the supersymmetric flavor and 
CP-problems by raising the masses of all scalars to several TeV. It was observed that for 
certain parameter combinations, the Higgs parameters have an infrared quasi-fixed point, 
thus making the weak scale soft parameter m# u of one of the Higgs doublets highly insensitive 
to the value of the other scalar masses. In this way, naturalness is preserved even for very 
large sfermion masses. 

Here, as a concrete example, the focus point scenario LCC2, as defined in the appendix, 
will be analyzed in detail. The collider phenomenology of this scenario has been studied 
previously in Ref. [27], and cosmological implications were discussed in Ref. [28]. In this 
report, the phenomeno logical discussion of the earlier works is extended, and the scenario 
compared to the baryogenesis case. 

3 Baryogenesis scenario 

The LHC experiments will be able to probe a light Higgs boson with Standard-Model-like 
couplings to the gauge bosons, as required by electroweak baryogenesis. The potential for 
discovery of light stops has been studied in detail for the LHC [14] and ILC [9, 15, 16]. 
If gluinos are not too heavy, the discovery of light stops in gluino decays via the process 
PP 99 — > ttt*t* : tttt has been shown to be effective for small stop-neutralino mass 
differences [14]. Stops can be discovered with more than five standard deviations through 
the decay of gluinos if the gluino mass is below 900 GeV, the stop-neutralino mass difference 
about 40 GeV, and all other squarks relatively heavy. It was also shown [14] that if the other 
squarks are not much heavier than 1 TeV, even in the stop-neutralino co-annihilation region 
with mi x — m^o < 30 GeV, stop discovery is possible in the decay of those squarks. However, 
further studies are needed to establish whether the stop-neutralino co-annihilation region 
can be explored at the LHC for first and second generation squarks with masses of several 
TeV. In this parameter region, the decay products of the stops are very soft and thus escape 
direct detection, so that only the final state products of the two top quarks from gluino 
decay are observable. Nevertheless, the signature of same-sign top quarks and large missing 
energy is likely sufficient for the discovery of a new physics signal, even if the stops are not 
directly identified. When discovered, the analysis of the kinematical decay distributions of 
the stops at the LHC can be used to extract information about the relationship between 

1 An analysis including two- loop corrections to the effective potential is in progress [26]. 
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stop, gluino and lightest neutralino masses [14], but an independent determination of the 
stop or neutralino mass seems difficult. 

At the 1LC, light stops can be discovered for stop-neutralino mass differences down to 
about 5 GeV, independent of other MSSM parameters [9]. In addition, precise measurements 
of the stop mass and mixing angle can be performed, with a typical stop mass error of about 
1%, assuming 80% beam polarization for electrons and 50% for positrons. 

3.1 Chargino and neutralino studies at the LHC 

In this kind of scenarios, the analysis of neutralinos and charginos is very difficult. Su- 
persymmetric particles are mainly produced through decay cascades of gluino and squarks, 
which can have large production cross-sections at the LHC. If however all squarks except 
the light stop are rather heavy, only the stop will be produced at sizeable rates and it will 
dominate the decays of the gluinos as well. However, a stop in accordance with electroweak 
baryogenesis is so light that it typically can only decay into the lightest neutralino, so that 
no information about the other neutralinos and charginos can be gained from stop processes. 

The only sfermion whose mass is not constrained by theoretical or experimental bounds in 
the MSSM baryogenesis scenario, is the right-chiral sbottom quark. The left-chiral sbottom 
is connected through SU(2) symmetry to the left-chiral stop, which needs to be very massive 
to satisfy the LEP Higgs bound [24]. The right-chiral sbottom, on the other hand, directly 
decays into a bino neutralino, with negligible branching ratios into other neutralinos or 
charginos. Note that the light sbottom state is expected to be almost completely right-chiral 
in this scenario, since large mixing between left- and right-chiral sbottom states is suppressed 
when the light mass eigenstate is supposed to be much lighter than the heavy state. 

Finally, charginos and neutralinos can also be produced directly through Drell-Yan-type 
processes. However, direct production processes of neutralino pairs have too small rates 
compared to the backgrounds at the LHC [29]. On the other hand, the production of a 
light chargino with a neutralino can have sizeable rates, but in a baryogenesis scenario, the 
chargino will mainly decay into the light stop. In the sample scenario BGEN (see appendix), 
the chargino branching fraction into stops and bottom quarks is more than 99.9%. As a 
result, the production of mixed neutralino-chargino pairs will lead to at most one lepton 
plus jets in the final state, instead of the typical tri-lepton signature for mSUGRA scenarios. 
Unfortunately, the one lepton plus jets signature is totally swamped by VU-boson background. 

Thus the only observable supersymmetric particles at the LHC are, depending on the 
region of parameter space, the gluino, the light stop, potentially the light sbottom, and the 
missing energy signature of the lightest neutralino. 

3.2 Chargino and neutralino studies at the ILC 

In the given scenario, stops, charginos and neutralinos can be studied precisely at a future 
linear collider. The prospective measurements for the stop are studied in detail in Ref. [9]. 
Here the chargino and neutralino phenomenology is investigated. The characteristic feature 
of the baryogenesis scenario is that while the charginos and neutralinos can be light, the first 
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Xi = xl xi xl 


X° = X? 

xl 

X3 

x2 


0.7 25.1 0.07 
0.4 62.0 0.005 


e+e" -> Xj 


xt = xt xi 


xj = Xi 

xl 


665 28 



Table 1: Tree-level production cross-sections in fb at ^/s = 500 GeV with unpolarized 
beams for the reference point BGEN. 

generation sleptons are heavy. A similar situation arises in focus point supersymmetry and 
has been studied in previous works [27,30]. Here the study of Ref. [30] is extended by includ- 
ing more observables, simulations for signal and background processes and by considering a 
CP-violating phase in the system. 

In general, future collider experiments will not provide enough independent measurements 
to extract all chargino and neutralino mass and mixing parameters in a model-independent 
way. Therefore any attempt to fully reconstruct a supersymmetric model from experimental 
data will rely on assuming some structure for that model. In particular, in the MSSM, the 
chargino and neutralino masses and couplings depend on five unknown parameters, namely 
Mi, M 2 , </> M and tan/3. In the MSSM, these parameters can be related to production 
cross-sections and masses. Thus with a sufficient amount of measurements for chargino and 
neutralino masses and cross-sections, all the underlying parameters can be extracted. 

For the sample scenario BGEN (see appendix), the chargino and neutralino masses at 
tree-level amount to 

m,o = 106.6 GeV, m,o = 170.8 GeV, m-± = 162.7 GeV, 

-M *2 q\ 

m,o = 231.2 GeV, m,o = 297.7 GeV, m-± = 296.2 GeV. 

At the ILC with y/s = 500 GeV, many neutralino and chargino states are accessible, see 
Tab. ^ As evident from the table, the most promising processes are e + e~ — > XiXii 
e + e~ — » X1X3 an d e + e~ — » x^xls- The production of mixed charginos pairs, e + e~ — ► x^X^j 
also has a sizeable cross-section, but is overwhelmed by background from e + e~ — > X1X1 ■ 

The relevant decay processes are summarized in Tab.|21 While the two lightest neutrali- 
nos dominantly decay into (virtual) Z bosons, the lightest chargino will decay into the light 
stop and a bottom. For light stops, the by far dominant decay mode is the loop induced 
process t — > cxI [31]. 

Chargino xt 

Production of light charginos, e + e~ — > XiXij nas a large cross-section of 665 fb at ^/s = 500 
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Sparticle 


Mass m [GeV] 


Width T [GeV] 


Decay modes 


X? 

xl 
xl 

xl 


106.6 
170.8 
231.2 

297.7 


0.00002 
0.11 

0.86 


xl^z*xl 100% 

xl^Zxi 98% 
^h Q xi 2% 

xl^zx\ 1% 
^Zxl «i% 
^Zxl 1% 
-> W± 94% 

->/i°xS 4% 
<1% 


xt 
xt 


162.7 
296.2 


0.24 
4.2 


X+ -> * 6 100% 
Xt - f6 84% 
-> w + Xi < 1% 
-> W + ^ 10% 
^Zx^ 6% 
^h°xt < 1% 



Table 2: Tree-level masses, widths and main branching ratios of the neutralino and chargino 
states at Born level for the reference point BGEN. 

GeV. It can be further increased by using beam polarization, P(e + )/P(e") = +50%/— 80%, 
where + stands for right-handed, — for left-handed polarization, yielding cr[e + e~ — > Xi~XT] = 
1760 fb. 

The dominant decay chain e + e~ — > xtxT ~* tt* bb — > cc bb XiXi leads to two charm 
jets and two bottom jets plus missing energy in the final state. The most important Stan- 
dard Model backgrounds arise from production of vector bosons in pairs (where the missing 
energy is generated by mismeasurements) and triples (where neutrino decays can lead to 
missing energy), as well as tt production. In addition, one needs to consider supersymmetric 
background from neutralino production. 

Signal and background events are simulated with the Monte-Carlo methods from Ref. [32], 
including full tree-level matrix elements and Breit-Wigner propagators for resonant in- 
termediate particles. The processes are generated on the parton level. Jet broadening 
through parton shower and detector effects are parameterized by smearing functions with 
lepton and jet energy uncertainty taken from [33]. Jets overlapping within a cone with 
AR = a/ (</>i — fa) 2 + (vi ~ V2) 2 < 0.3 are combined into one jet, where 0j and rji are the 
azimuthal angle and rapidity of jet %. Similarly, a lepton lying within a jet is combined into 
the jet. Leptons and jets outside the central region of the detector have a higher likelihood 
of mistag and get inflicted by large two-photon background. Therefore leptons within an 
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angle of | cos 6\ < 0.95 around the beam line and jets with | cos 6\ < 0.90 are discarded. After 
these numerical adjustments the remaining isolated jets and leptons define the signature of 
the simulated event. 

Background from two-photon interactions was not simulated for this study, since they 
typically lead to very soft hadronic events, and can be removed by a cut on the total trans- 
verse momentum, pl ot > 12 GeV [9]. 

The other backgrounds are reduced by the following simple kinematic cuts: 

• Each event must contain four hadronic jets and no isolated lepton and pl ot > 12 GeV 
(see above). 

• For each combination of two jets, the invariant mass is required to be different from the 
Z mass, \m,jj — m z | > 10 GeV. This strongly reduces background from gauge bosons 
and neutralinos. 

• Since backgrounds from gauge bosons tend to increase in the forward and backward 
detector regions, they can be further reduced by a cut on total momentum angle, 

I COs6 WJ = |Plong,tot/Ptot| < 0.9. 

• All Standard Model backgrounds are cut down by requiring large missing energy, Ifi > 
100 GeV. 

• Since two bottom jets are expected in the signal, the signal-to-background ratio is 
improved through b-tagging. Following [34] , it is assumed that the b tagging efficiency 
is 90%, with a mistag probability of light flavors of 10%. B-tagging works on all 
backgrounds except tt. For the other backgrounds the four jets mainly originate from 
W decays, so that two of the jets could be charmed. Charm jets have a higher b- 
mistagging probability than light flavors, but since the signal also always contains two 
charm jets, the resulting tagging power between signal and gauge boson backgrounds 
is governed by the discrimination between bottom and light flavors. 

With this selection procedure, the remaining background is very small, about 3 fb, whereas 
the total signal efficiency is 18%. Including an overall systematic acceptance of 90%, the 
resulting signal efficiency is 16%. With a total luminosity of 250 fb -1 the statistical error 
for the cross-section measurements is 

5a RL [ X fxl] = 0.37%, 6<j LR [xtxl] = 1-6%, (2) 

where RL and LR stand for the polarization combinations P(e + )/P(e _ ) = +50%/— 80% and 
-50%/+80%, respectively. 

The spectrum of the decay products can be used for a determination of the chargino 
and neutralino mass. The spectrum of the invariant mass of the b and c jet from the decay 
xl — > tb — > cbxi has a characteristic upper edge at 

m bc,max = ™bc,max = ™-± ~ m^O. (3) 
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Figure 1: Distribution of the bottom-charm invariant mass m^ c from chargino decay and 
a simple fit to the upper edge. 

However, with a four jet final state the problem remains to identify the pairs of two jets each 
that belong to the decay of one chargino. The bottom jets can be identified as the jets with 
the highest b tagging likelihood. However this still leaves a twofold ambiguity to combine 
a b jet with a c jet. Here the following strategy is adopted: since a pair of jets originating 
from different charginos tend to have a larger invariant mass than a pair originating from 
the same chargino, always the pair with the lower invariant mass is selected for the mass 
measurement. 

The resulting invariant mass distribution in shown in Fig. A very crude fit of the 
upper threshold gives 



stemming from the jet energy scale. Based on studies at LEP for W boson pair production 
[35], the jet energy scale error is assumed to be 0.4%. 

Since this measurement would only yield the difference between chargino and neutralino 
mass, an independent measurement is necessary for the determination of absolute values of 
both masses. A very precise direct determination of the chargino mass can be obtained from 
a threshold scan. By measuring the chargino cross-section at a few center-of-mass energies 
near the pair production threshold, the onset of the pair production excitation curve can 
be reconstructed and used for a mass determination. By including two points below the 
threshold, the background can be analyzed and extrapolated in a model-independent way. 
The result of the analysis for five scan points with 10 fb _1 luminosity each, using the same 
methods and selection cuts as above, is depicted in Fig. The result of a fit to the threshold 
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Figure 2: Threshold scan for chargino pair production using five scan points with 10 fb^ 1 
each. 

excitation curve is 

m-± = | [Vs\ thi = (163.02 ± 0.04 ± 0.04) GeV, (5) 

where the first error is statistical and the second is the systematic error from the beam 
energy calibration. Combining eqs. (j3J) and (J3J), the masses are obtained as follows, 

m_± = (163.02 ± 0.06) GeV, m^o = (106.1 ± 0.3). (6) 
Neutralinos x® an d X° 

Neutralino production can have sizeable cross-section for mixed pair production of X1X3 
and X2X3- The signal for X1X3 is characterized by two leptons or two jets plus missing 
energy, which receives very large backgrounds from W + W~ and ZZ production. Even after 
reducing the backgrounds with appropriate cuts, the remaining level is still problematic for 
any precision measurement. Therefore, in the following only X2X3 is considered. 

This process has a sizeable cross-section, a[e + e~ — > X2X3] = 62 fb for unpolarized beams. 
Using right-handed polarization for the e~ and left-handed polarization for the e + beam, 
P(e + )/P(e~) = -50%/+80%, the signal is enhanced to 68.9 fb, while all Standard Model 
backgrounds, which are dominated by left-chiral SU(2) interactions, are reduced for this 
beam polarization combination. 

The X3 almost always decays into a Z boson and lightest neutralino Xii while the 
due to the small mass difference m^o — m^o decays through a virtual Z into two quarks or 
two leptons. In order to improve the statistical significance of the signal, the hadronic decay 
modes of the Z are considered, since the branching ratio of Z to leptons is very small. Then 
the final state is characterized by four (light quark) jets and missing energy. Standard Model 
backgrounds arise from processes with two and three vector bosons and from ti production. 
In addition, chargino production as discussed above, constitutes another difficult background. 
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Signal and background are simulated and analyzed as detailed above. The following cuts 
have been applied to select the signal out of the background: 



• Similar to the chargino analysis, a four jets signature and the same cuts on cos# Ptot 
and Ip are used. 

• The invariant mass of two jets has to be equal to the Z boson mass, \m,j a j b — toz| < 10 
GeV, whereas the invariant mass of the other two jets has to be smaller than toz, 
mz — rrij c j d > 10 GeV. All combinatoric pairings of the four jets are tried for this 
purpose, and if one combination meets these requirement, the event is kept. 

• The Background from charginos generates rather soft jets due to the small stop mass. 
Thus a cut on the transverse momentum, pj: ot > 50 GeV is effective against this 
background. 

• Finally the chargino background is further reduced by a b-quark veto. As a conse- 
quence, here the signal is restricted to light quark decay channels of the Z bosons 
only. 

After this selection procedure, the remaining background is around 0.2 fb, but a good signal 
efficiency of 26% is achieved. Including a general systematic acceptance of 90%, the resulting 
signal efficiency is 24%. With a total luminosity of 250 fb -1 for the polarization combination 
P(e + )/P(e~) = — 50%/+80%, the statistical error for the cross-section measurement is 

SvrlIxIxI] = 1-6%. (7) 

For the opposite polarization combination P(e + )/P(e~) = — 50%/+80%, the precision of the 
neutralino cross-section measurement is much worse due to larger backgrounds, and does not 
improve the global MSSM analysis. 

Information about the neutralino masses can be extracted from the decay product distri- 
butions. In the Xg — > Zxl decay, the energy spectrum of the Z is rather flat, with distinct 
lower and upper endpoints at 

£min, m ax,3 = (™*° " + ™|°™Z " 4™Z + m \f + ™\ S 

xg vsv " " (8) 



2/2 2 

m~o (to -o - to -o + s 

Al x A3 A2 



) T yj K m \o, m 2 -o, to|) A(to|o, m 2 ~o, s)) , 



with A (a, b, c) = a 2 + b 2 + c 2 — lab — lac — 26c. The energy of the Z boson can be deduced 
from the energy of the jet pair that combines to the Z invariant mass. 

For the \\ decay, the kinematics are different since only a virtual Z is involved in the 
process. Still the energy spectrum of the two final state jets has a distinct upper endpoint 
given by 

to-o — to-o — 2m 2 ov / i + s 
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Figure 3: (a) Energy distribution for jet 
pair in \\ Zxl decay, (b) energy distribu- 
tion for jet pair in —* Z*Xi decay, and (c) 
invariant mass of the latter, including simple 
fitting curves. 



In addition, the invariant mass of the same two jets is bounded from above by the mass 
difference of the \2 an d X?, 

"ijj,max,2 = rn^o - m^o. (10) 

The simulated distributions after cuts and fits to the endpoints are summarized in Fig. 01 
From a combination of the different fits, and including the jet energy scale error as for the 
charginos, one obtains 

5m^ = 1.3 GeV, 5m^o = 1.5 GeV, 5m^o w 3.3 GeV. (11) 

By feeding in the more precise determination of the lightest neutralino mass in eq. (|B]L these 
numbers are improved to 

5m^o = 0.6 GeV, 5m^o w 2.0 GeV. (12) 

The mass measurements are not sufficient to extract the fundamental supersymmetry 
parameters without ambiguity. Therefore the additional information from measurements of 
the chargino and neutralino cross- sect ions are important. The MSSM model parameters can 
be extracted from a combined fit to all observables. 

As a result of very heavy selectrons and sneutrinos, the t-channel contribution to the 
production processes are essentially switched off. However, the non-observation of these 
particles at colliders puts only a relatively mild lower bound on their masses of about 500 
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GeV. Here the values of the sleptons masses are kept arbritrary, and instead indirect bounds 
on the masses are derived from the fit (see also Ref. [30]). 

The chargino and neutralino cross- section and mass measurements are combined in a x 2 
fit, thus taking into account all correlations in the underlying parameters that enter in the 
observables. The fit results with one standard deviation errors are 

Mi = (118.8 ± 0.4) GeV, 
M 2 = (225.0 ± 0.9) GeV, 
\n\ = (225.0 ± 1.2) GeV, 
|^| < 0.6, 

The resulting constraint on the phase (p^ is rather poor, since none of the included observables 
is directly CP-sensitive, and M is strongly correlated with other parameters, in particular 
tan f3. 



tan j3 



c+0.5 
J -0.7> 



m,~, > 5 TeV, 



(13) 



Wej, > 1.5 TeV. 



CP violation 

The effect of CP violation can be uniquely studied in observables that are directly CP- 
sensitive. In the chargino and neutralino sector, triple products of kinematic momenta have 
been identified as useful for that purpose [36-38] . 

In this work, hadronic decays of neutralinos and charginos are investigated, 

e + e" -> x S iXj s s = 0, ±. 

^ JaJbXl 

From the momenta p^b of the two jets one can construct the T-odd triple product 

T = p c -(p^x p h ). (15) 

and the T-odd asymmetry 

a[T>Q}-a[T<0] 
' A ~ a[T>0] + a[T<0Y 1 bj 

The experimental investigation of this triple product requires the identification of the charge 
of the jets. While this is not possible for light-quark jets on a event-by-event basis, it can 
be achieved for larger samples of jet events on a statistical basis, see e.g. Ref. [39]. 

Using the same Monte-Carlo techniques as explained above, the asymmetry is calcu- 
lated both for chargino and neutralino production. In the given scenario, it is found that 
^[XiXil = 0.3% and -A^xl] = 0.9%. Thus in both cases the expected magnitude is of the 
same order as the statistical error of the cross-section measurements, see eqs. ©, ©• 

For larger values of M , larger triple product asymmetries can be obtained, however 
such models are already highly constrained by electric dipole moment limits. Therefore a 
measurement of the CP-violating phase in the fi parameter appears to be hopeless at the 
ILC. These findings are in agreement with previous theoretical studies [37]. 
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3.3 Cosmological implications 

The discovery of light scalar top quarks, in conjunction with a Standard-Model-like Higgs 
boson with a mass near 120 GeV, would be a strong indication that electroweak baryoge- 
nesis is the mechanism for the generation of the baryon asymmetry. At the same time, 
supersymmetry could also explain the existence of dark matter in the universe, based on the 
co-annihilation mechanism. In order to confirm this exciting picture, the relevant supersym- 
metry parameters have to be measured accurately. 

One needs to (i) determine that the light stop is mainly right-chiral to contribute ap- 
propriately to the electroweak phase transition while being in agreement with electroweak 
precision measurements, (ii) check that the masses and compositions of the gauge/Higgs 
superfield sector are compatible with the values required for the generation of the baryon 
asymmetry, and (iii) compute the dark matter annihilation cross-sections and the relic abun- 
dance so as to compare with cosmological observations. If stop-neutralino co-annihilation is 
relevant it is important to determine the stop-neutralino mass difference very precisely. 

Baryogenesis 

In the given scenario, the first information about the necessary constituents for baryogenesis 
in the MSSM can be obtained at the LHC, namely the discovery of a light Higgs boson 
and possibly evidence for a light stop. At the ILC, the picture can be rendered much more 
precisely by determining the mass and mixing angle of the light stop accurately. 

Using the computations of Refs. [5,20] for the electroweak phase transition, it is found 
that the experimental uncertainty of the stop parameters allows to determine the strength 
of the phase transition to better than 10%, 



Note however that this is only a parametric error, and there are also theoretical uncertainties 
involved in the computation. 

Besides the strength of the electroweak phase transition, the second crucial ingredient for 
electroweak baryogenesis is the generation of the baryon asymmetry through CP-violating 
processes. As explained in section |21 in the MSSM these CP-violating contributions can 
be described by a complex phase of the parameter /i. However, the findings of the previous 
sections show that for typical values of the CP- violating effects in neutralino and chargino 
observables at the ILC are too small to be observed. It would only be possible to derive an 
upper bound of |0 M | < 0.7 at the 90% confidence level, which would leave the question of 
the baryon asymmetry still undecided. 

Dark matter 

As discussed widely in the literature (see e.g. [9,28,40]), collider experiment data can be 
used to compute the expected cosmological dark matter relic density. In the scenario studied 
here, the relic density is governed by the neutralino-stop co-annihilation mechanism, and thus 
crucially depends on the stop-neutralino mass difference. It was found [41] that this mass 
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Figure 4: Computation of dark matter relic abundance JIcdm^ 2 for the scenario BGEN, 
taking into account estimated experimental errors for stop, chargino, neutralino and Higgs 
sector measurements at future colliders. The dots correspond to a scan over the la (Ax 2 < 1 ) 
region allowed by the experimental errors, as a function of the measured stop mass, for a 
mass measurement error of 1.2 GeV (gray dots) and 0.3 GeV (black dots). The original 
scenario used as input is indicated by the red (light gray) star. The horizontal shaded bands 
show the la and 2a constraints on the relic density measured by WMAP and SDSS. 

difference cannot be extracted with good precision directly from stop decay distributions, 
because of radiation effects and limited statistics, but the best determination is derived 
from independent mass measurements of the stop and neutralino. While the stop mass 
and mixings can be extracted from precision measurements of the stop cross-section [9] , the 
neutralino properties have to be determined from independent observables as discussed in 
this report. 

The relic dark matter density is computed with the codes described in Ref. [12,22]. 
The analysis is based on observables in the stop and neutralino/chargino sectors. For the 
neutralinos and charginos, the estimated errors taken from the previous chapter, see eqs. (J2J), 
©, ©, jUJ)- The stop measurement errors are derived from the study of Ref. [9] with the 
result: Sm^ = 1.2 GeV and | cos0 t ~| < 0.077. 

In total, the relevant parameters used as input are m^, cos 6^-, Mi, M 2 , M , tan/?. 
It has been checked that the dependence on Higgs and slepton parameters is negligible for 
this scenario. The mass of the heavier stop t 2 is too large to be measured directly, but it is 
supposed that a limit of mi 2 > 1000 GeV can be set from LHC searches. 

The expected experimental errors are propagated and parametric correlations are taken 
into account by means of a \ 2 scan. However, theoretical errors due to missing radiative 
corrections are not taken into account in the calculation of the dark matter density. Fig. 0] 
shows the result of a scan over the MSSM parameter space for the scenario BGEN. The 
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scattered gray dots indicate the region allowed by the collider experimental uncertainty, as 
a function of the measured stop mass. The range of the horizontal axis is constrained by the 
error in the stop mass measurement, m~ tl = (122.5 ± 1.2) GeV. The horizontal bands depict 
the relic density as measured by WMAP and SDSS [1,2] with one and two standard deviation 
errors. At la level, the astrophysical observations lead to 0.103 < QcBMh 2 < 0.116. 

In total, using the collider measurement simulations, the relic density can be predicted 
to 0.082 < ^cdm^ 2 < 0.139 at the la level. Thus the overall precision is of the same 
magnitude as, but worse by roughly a factor 4 than the direct WMAP/SDSS determination. 
The uncertainty in the theoretical determination is dominated by the measurement of the 
ti mass with an error of 1.2 GeV. This shows that even in a scenario with heavy sfermions, 
which is unfavorable for future collider experiments, the achievable precision in the analysis 
of the chargino and neutralino sector is very high and not a limiting factor for cosmological 
interpretations. 

First results of an optimized threshold scan method indicate that the precision for m~ tl 
can be improved to about 0.3 GeV [42]. The advantage of the threshold scan method is 
the small influence of systematic errors, since it makes use of the shape of the cross-section 
as a function of the center-of-mass energy, instead of absolute cross-section measurements. 
However it is limited by small statistics near the threshold. The best accuracy can be 
achieved by combining a measurement near the threshold, where the cross-section is most 
sensitive to the stop mass, with a measurement at higher energies, where the cross-section is 
larger [42]. With a stop mass error of Sm^ = 0.3 GeV, the relic density could be computed 
much more precisely, yielding the result 0.099 < ^cdm^ 2 < 0.122 in the scenario BGEN. 
This precision is very comparable to the direct WMAP/SDSS determination, as indicated 
by the black dots in Fig. HJ 

4 Focus-point scenario 

4.1 Chargino and neutralino studies at the LHC 

In a focus-point scenario, all sfermions are very heavy, but gluinos g can be sufficiently light 
to be generated with large cross-sections. The gluino decays can proceed via some neutralinos 
and charginos, giving a clear and distinct leptonic discovery signature [29]. For example, 
a signature with same-sign leptons and missing transverse energy can be interpreted as a 
decay cascade of the gluinos via charginos [43]. Besides, the decay \2 ~~ * Z*Xi ~ > l + l~Xi 
allows to determine the mass difference m 21 = m^o — m^o from the invariant mass 
distribution with a precision of about 6m 2 i ~ (0.5 GeV) x exp(— nig/ (200 GeV) + 3.5), 
where the exponential factor with nig approximately accounts for the decrease of the gluino 
cross-section with larger gluino masses. 

4.2 Chargino and neutralino studies at the ILC 

In the scenario LCC2, among the superpartners, only charginos and neutralinos can be 
studied at a future linear collider, whereas all sfermions are too massive to be accessible 
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~o ~ n ~- n ~ o 

Xi = X2 X3 X4 


X° = X? 

xl 

X3 

x2 


8.0 46.3 0.01 
8.7 107.2 0.07 
6.2 22.7 


e+e" -> Xj 


xt = xt xf 


xj = Xi 

xl 


534 23 



Table 3: Tree-level production cross-sections in fb at ^/s = 500 GeV with unpolarized 
beams for the reference point LCC2. 

either at the LHC or ILC This situation arises always in focus point supersymmetry and 
has been studied in previous works [27,30]. Here the study of Ref. [30] is extended by 
including different observables and simulations for signal and background processes. 

As in the previous section, the general procedure for reconstructing the underlying MSSM 
parameters relies on the assumption of the MSSM structure for the chargino and neutralino 
mass matrices. The parameters Mi, M 2 , M and tan/3 are extracted from a fit to mass 
and cross-section measurements. 

For the scenario LCC2 (see appendix), the chargino and neutralino masses at tree- level 
amount to 

m,o = 107.7 GeV, m,o = 166.3 GeV, m~± = 159.4 GeV, 

*1 *2 Xl (IS) 

m,o = 190.0 GeV, m,o = 294.4 GeV, m-± = 286.8 GeV. 

As pointed out above, the LHC has a good possibility to find evidence for and measure 
some properties of neutralinos and charginos. At the ILC with ^Js = 500 GeV, many 
neutralino and chargino states are accessible, see Tab. El The most promising processes are 
e + e~ — > Xi~X;i~5 e + e~ — > X1X3 an d e + e~ — > X2X3; similar to the baryogenesis scenario. 

Tab. |U gives the relevant decay modes and branching fractions. In this scenario, the 
mass differences between the lighter neutralino and chargino states are rather small, thus 
allowing only decays of the light states through virtual gauge bosons. 

Chargino xt 

Production of light charginos, e + e~ — > X1X1 > nas a large cross-section of 534 fb at yfs = 500 
GeV, which can be further increased to 1345 fb by using beam polarization, P(e + )/P(e~) 
= +50%/-80%. 

In the focus-point scenario, charginos almost always decay through virtual W bosons, 
e + e~ — > X1X1 ~^ W + *W~* X1X1 ~^ \ ±v i <?<?' Xi Xi where it is useful to consider a final 
state where one W* decays hadronically into quarks, while the other decays into a lepton 
I = e, fi and a neutrino. This final state allows a precise kinematical measurement, since the 
two intermediate virtual W bosons can be disentangled, as opposed to the purely hadronic 
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Sparticle 



X? 

xl 



xl 



xt 



Mass m [GeV] 



107.7 
166.3 
190.0 

294.4 



159.4 
286.8 



Width T [GeV] 



0.000016 
0.0013 

0.81 



0.0002 
0.73 



Decay modes 



Xil 

xl 



xl 



Z*X°i 

z*xi 
zx\ 

zxl 
zx% 

W ± xf 
h°x°2 



Xi 

xt 



. w + xi 

■W + x°2 

. w + xi 

■ZxX 

h°xt 



100% 
100% 

< 1% 

< 1% 
1% 
1% 

94% 

< 1% 
3% 



100% 

3% 
42% 
11% 
36% 



Table 4: Tree-level masses, widths and main branching ratios of the neutralino and chargino 
states at Born level for the reference point LCC2. 



final state. As before, Standard Model background from vector bosons and tt, as well as 
super symmetric background from neutralino production is taken into account. 

Signal and background events are simulated with the same Monte-Carlo methods as 
elaborated above. The following selection cuts are applied to reduce the backgrounds: 

• Each event must contain two hadronic jets, one isolated lepton and pl ot > 12 GeV. 

• The invariant mass of the two jets is required to be smaller than the W mass, m-yv — 
rrijj < 10 GeV, to reduce Standard Model background. 

• All Standard Model backgrounds are cut down by requiring large missing energy, Ifi > 
100 GeV, and the direction of the missing momentum to be in the visible detector 
region, I cos#„ . I < 0.8. 

O 7 I /'miss I 

• Since a large background comes from W + W~ production, the invariant mass of the 
lepton and the missing momentum, mip miss tends to increase near the W mass. Thus 
by requiring mi Piniss > 150 GeV, that background is reduced. 

After these cuts, the remaining background is very small, but 39% of the signal is retained. 
With an overall systematic acceptance of 90% and a total luminosity of 250 fb _1 the statistical 



18 



error for the cross-section measurements is 

Sa RL [ X tXi] = 0.6%, 5a LR [xtxi] = 1-7%, (19) 

where RL and LR stand for the polarization combinations P(e + )/P(e~) = +50%/— 80% and 
-50%/+80%, respectively. 

The xt an d Xi masses can be constrained from distributions of the decay products. The 
summed energy of the two jets from the decay Xi ~^ W^Xi — * QQ'Xi have a well-defined 
upper endpoint at 

Moreover, the invariant mass of the same two jets is bounded from above by 

"ijj.max = m^± - m^o. (21) 

Simple fits to the upper tails of the energy and invariant mass distributions give 

Pjj.max = (120.1 ± 0.6 ± 0.5) GeV, m jjimax = (51.8 ± 0.14 ± 0.2) GeV, (22) 

compared to the model input values -Ejj.max = 120.2 GeV and m~± — m^o = 51.7 GeV. Here 
the first error is the statistical uncertainty of the fit, whereas the second error accounts for 
the 0.4% jet energy scale uncertainty. These measurement translate into the following results 
for the masses: 

m_± = (159.4 ± 1.0) GeV, m^o = (107.7 ± 0.9) GeV, (23) 

with a large correlation between the two masses. With an independent determination of 
one of the masses from a different observable, both mass determinations can be improved 
substantially. 

As before, a threshold scan of the chargino cross-section yields a much more precise 
determination of the chargino mass. Using five scan points with 10 fb _1 luminosity each and 
the same selection cuts as above, a fit to the threshold excitation curve gives 

m xf = \ [v^thr = (159.38 ± 0.04 ± 0.04) GeV, (24) 

where the first error is statistical and the second is the systematic error from the beam 
energy calibration. Combining eqs. (|22jl and (|24jl . the masses are obtained as follows, 

m _± = (159.38 ± 0.06) GeV, m^o = (107.7 ± 0.21). (25) 
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Neutralinos x® an d X3 

As for the baryogenesis scenario, only the neutralino production process X2X3 yields a good 
signal-to-background ratio. This process has a sizeable cross-section, a[e + e~ — > X2X3] — 
107 fb for unpolarized beams. Using right-handed polarization for the e~ and left-handed 
polarization for the e + beam, P(e + )/P(e~) = — 50%/+80%, the signal is enhanced to 129 fb, 
while all Standard Model backgrounds are reduced for this beam polarization combination. 

Both the X2 an d X3 almost always decay through a virtual Z boson to a pair of leptons 
or quarks and the lightest neutralino. As a consequence, the decay products of the two 
neutralinos cannot be distinguished by kinematical constraints. Therefore, here the semi- 
leptonic decay channel is investigated, with one Z decaying leptonically, and the other Z 
decaying hadronically. Thus the final state is characterized by two charged leptons, two jets 
and missing energy. As for the baryogenesis case, the relevant Standard Model backgrounds 
arise from processes with two and three vector bosons and from ti production. 

The selection cuts to extract the signal from the background are very similar to the 
baryogenesis case: 

• Each event must contain two hadronic jets and two isolated lepton and pl ot > 12 GeV. 

• The invariant mass of either the two jets or the two leptons has to be smaller than mz, 
rriz — rrijj > 10 GeV and mz — mu > 10 GeV. 

• As in section E2J it is required that \cos9 ptot \ < 0.9 and Ifi > 100 GeV and a b-quark 
veto again ti background is applied. 

A resulting signal efficiency of 47% is obtained. With a total luminosity of 250 fb -1 for the 
polarization combination P(e + )/P(e _ ) = — 50%/+80%, the statistical error for the cross- 
section measurement is 

SvRtixlxi] = 2-7%. (26) 

For measurements of kinematic distributions, it is useful to restrict oneself to the two leptons 
in the final state only, which can be measured more cleanly and precisely than jets. The 
summed energy of the lepton pair has a distinct upper endpoint, depending from which of 
the two neutralinos they originate, 

771-0 — m 2 ~o — 2m 2 - y/s + s 

771 X2 X3 Xl (O f 7\ 
-&ll,max,2 — J= , {*() 

777-0 — 777 - — 2 771-0^/5 + S 
Al,max,3 — • {^O 

In addition, the invariant mass spectrum of the same two leptons has an endpoint at the 
mass difference between the intermediate and the final state neutralino, 

W U) max,2 = m^O ~ TTl^o, m U)maXi3 = 777^0 - 777^0 (29) 

The simulated distributions after cuts and fits to the endpoints are summarized in Fig. El For 
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Figure 5: (a) Energy distribution for lepton pair in X23 ^Xi decays, (b) invariant mass 
of the latter, including simple fitting curves. 



the masses in Tab.EJ the upper endpoint of the E\\^ spectrum is smaller than the maximum 
of the En >3 spectrum originating from the As a result, the endpoint £ii !max ,2 gets diluted 
in the continuum of the En )3 distribution. Therefore, here only the endpoint E^ max 3 is 
fitted, but no numerical value for En jmaXj 2 is obtained. The m\\ shows two characteristic 
edges corresponding to the contribution from the X2 an d the 

Combining fits to these three endpoints, and including a lepton energy scale error of 0.1%, 
one obtains 



Al 



2.9 GeV, 



X.2 



-tH GeV, 



A3 



-+t 8 GeV. 



(30) 

By feeding in the more precise determination of the lightest neutralino mass in eq. (j'25jl . 
these numbers are improved to 



Sm^o =±™ GeV, 



Sm^o =t 2 2 6 3 GeV. 



(31) 



From a fit to the cross-sections, not only constraints on the chargino and neutralino 
parameters are obtained, but also indirect bounds are derived for the masses of the heavy 
sleptons appearing in the t-channel production contribution (see also Ref. [30]). 

The result of the combined fit, taking into account chargino and neutralino cross-section 
and mass measurements, is 



Mi = (123.1 
M 2 = (237.6 
= (178.6 
\d>J < 0.6, 



+0.4\ 
-0.3J 
0.7^ 
l.lJ 
0.5\ 
0.5J 



GeV, 
GeV, 
GeV, 



tan/3 = 10; 



h0.8 
-1.4) 



2.6 TeV < m c 



m 



< 4.8 TeV, 
> 1.2 TeV. 



(32) 



where the errors indicate one standard deviation (la) uncertainties. The precision for the 
underlying MSSM parameters is comparable to the findings of Ref. [30], although this com- 
parison is precarious since the focus-point scenario studied there is different from the point 
LCC2. The underlying scenario LCC2 does not contain any CP-violating phases, but nev- 
ertheless a non-zero phase of the /i parameter was allowed in the fit. It turns out that the 
mass and cross- section observables give only a rather poor constraint on M . 
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Figure 6: Computation of dark matter relic abundance S^cdm^ 2 for the scenario LCC2, as 
a function of the lightest neutralino mass m^o. The black dots are the la allowed region. 
The original scenario used as input is indicated by the red (light gray) star. The horizontal 
shaded bands show the la and 2a constraints on the relic density measured by WMAP and 
SDSS. 




4.3 Cosmological implications: dark matter 

In focus-point supersymmetry, the annihilation cross-section is driven by a sizeable higgsino 
component of the lightest neutralino. As a result, the neutralinos can annihilate efficiently 
into electroweak gauge bosons. In order to scrutinize this scenario, the neutralino properties 
need to be determined with high precision. For the sfermions, only indirect bounds can be 
obtained, but this is sufficient to establish the fact that they effectively decouple from the 
annihilation cross-section. 

The relic dark matter density is computed with DarkSUSY 4.1 [44]. Including the 
expected precision for the hypothetical neutralino and chargino measurements discussed 
in the previous section, see eqs. (fTUjl. (|23j) . (J2S])> (GUJ), the accuracy for the dark matter 
determination is obtained from a x 2 scan. It is assumed that a limit of irtq > 1000 GeV for 
the squark masses can be set from non-observation of these particles at the LHC While in 
many scenarios it is possible to explore much larger squark masses, a bound of 1000 GeV 
is sufficient to guarantee that the effect of the squarks on the dark matter annihilation is 
negligible. 

Fig. depicts the result of the scan for the scenario LCC2. The scattered black dots 
indicate the region allowed by the collider experimental uncertainty, as a function of m^o. 
As pointed out above, although the LCC2 scenario is CP-conserving, a CP- violating phase for 
/i was allowed in the fit to the mass and cross-section measurements. However it was checked 
that when constraining M to zero in the fit, the results remain essentially unchanged. The 
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reason for this is that the annihilation cross-section only mildly depends on the CP-phase, 
but much more strongly on the neutralino mass parameters, which are directly measured. 
The precision of the predicted relic density f2cDM^ 2 is a very remarkable 2.5% at la level. 

The computed relic density depends sensitively on the lightest neutralino mass ffi^o. 
Compared to Refs. [27,28], the neutralino mass determination is improved in this analysis 
by the inclusion of the chargino threshold scan. As a result, the dark matter computation 
from the hypothetical collider measurements is also more accurate, with 5[Qcr>Mh 2 } ~ 2.5% 
instead of 5[n CD Mh 2 } ~ 8% in Ref. [28]. 



5 Impact of positron polarization on ILC results 

Throughout the previous sections, is was assumed that both the e~ and the e + beams at 
the ILC can be polarized, with polarization degrees of |P(e + )| = 50% and |P(e~)| = 80%. 
However, since the realization of positron polarization poses a serious challenge for the 
accelerator design, it is interesting to investigate how much the results presented above 
depend on it. As an extreme case, in this section the situation of zero positron polarization is 
studied, while as before 80% polarization is assumed for the electron beam, \P(e + ) \/\P(e~) | = 
0%/80%. The availability of e~ polarization ensures that all observables that are considered 
in the previous sections for the ILC still exist in this case. In particular, supersymmetric 
production cross-sections can be measured for two different polarization values, P(e~) = 
±80%. However, a loss of accuracy can result from the absence of positron polarization. 

In the baryogenesis scenario BGEN, the uncertainty of the determination of the elec- 
troweak phase transition strength and of the dark matter density is dominated by the error 
in the stop mass. The stop mass can be extracted from measuring the stop production 
cross-section for two different beam polarization combinations. Including 50% positron po- 
larization [9], an error of Sm^ = 1.2 GeV is found with this method, while without positron 
polarization, the error is about 20% larger, Sm^ = 1.4 GeV. Similarly, the accuracy for the 
chargino and neutralino parameters M 1? M 2 and [i is reduced by roughly 20% with respect 
to the values in eq. (fT3~j) when positron polarization is absent, but these parameters have a 
smaller impact on cosmological quantities. As a result of the larger stop mass error, a bigger 
uncertainty for the prediction of the dark matter density is obtained, f2cDM^ 2 = O.IOO^q;^ 2 , 
compared to loiol? with 50% e + polarization. As mentioned at the end of section l3~31 the stop 
mass can also be determined more precisely from a second cross-section measurement near the 
stop pair threshold. This method is less sensitive to the beam polarization, and with zero e + 
polarization the errors are only slightly larger: Sm^ = 0.32 GeV and ^cdm^ 2 = 0.109^q;^25 
instead of 5m h = 0.30 GeV and fi C DM^ 2 = 0.109^;^. 

For the focus-point scenario LCC2, the low-energy phenomenology is governed by the 
chargino and neutralino states. With 0% instead of 50% position polarization, the uncer- 
tainty of the parameters in eq. (|3^j) is larger by roughly 20% to 30%. The upper bound that 
can be extracted for the sneutrino mass from the chargino cross-section measurements gets 
much weaker, m/> c < 13 TeV, since the derivation of this parameters is very delicate and 
requires high precision. The larger errors in the chargino and neutralino parameters also 
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translate into a larger error on the predicted dark matter density of 3.1%, compared to 2.5% 
with 50% position polarization. Nevertheless, the possibility to compute the dark matter 
density with about 3% precision from collider data is still a very impressive result. 

6 Conclusions 

In various theoretical supersymmetry frameworks, the mass of most or all scalar fermions 
is constrained to be very large, of the order of at least a few TeV. Typical examples are 
electroweak baryogenesis in the MSSM, focus-point supersymmetry and split supersymmetry. 
Focusing on two benchmark points, one electroweak baryogenesis scenario and one focus- 
point scenario, it was studied how such a supersymmetry scenario can be explored at future 
colliders. It was found that the LHC can make discoveries of several superpartner particles 
in the focus-point case, while in the baryogenesis scenario a new physics signal is likely to be 
observable at LHC, but the identification of the contributing superpartner particles is very 
challenging. In both scenarios, the ILC can perform precision measurements to determine 
the properties of these particles and set bounds on other unobserved states. 

The analysis was performed including realistic signal and background computations and 
simple evaluations of statistical and systematic experimental errors for mass and cross-section 
measurements. The study was based on the full MSSM, i.e. without assuming a specific 
mechanism or pattern for supersymmetry breaking parameters. 

Using these results, the cosmological implications for electroweak baryogenesis with light 
stops and stop-neutralino co-annihilation were investigated. It turns out that the collider 
data helps to elucidate the strength of the electroweak phase transition, while the manifesta- 
tion of the CP- violating source responsible for the baryon asymmetry remains unconstrained. 
Furthermore, the collider measurements can be used to compute the relic dark matter density. 
By determining the stop and lightest neutralino masses, the stop-neutralino co-annihilation 
process can be strongly constrained and the dark matter density predicted with a precision 
of the same order as current astrophysical results. Refinements in the determination of the 
stop mass can improve this result significantly. 

In the focus-point scenario, it was found that very precise measurements of the accessible 
neutralino and chargino states can be performed at the ILC, combining mass measurements 
from distributions and threshold scans together with cross-section measurements. This allows 
to set constraints on the slepton mass scale and to compute very accurate predictions for the 
relic density at the per-cent level. Similar conclusions also apply for split supersymmetry if 
the lightest chargino and neutralino particles have masses of a few 100 GeV. 

In both scenarios, the dependence on positron beam polarization was investigated. As- 
suming 80% e~ polarization, and comparing the case of 50% e + polarization to zero e + 
polarization, it is found that all elements of the analysis can be performed similarly, but 
roughly 20-30% precision in the relevant supersymmetry parameters and derived cosmolog- 
ical quantities is lost without e + polarization. 

The present work has been performed using tree-level formulae and cross-sections for the 
ILC analysis and the computation of the dark matter annihilation rate. For the expected 
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experimental precision, however, radiative corrections are important and introduce a depen- 
dence on other supersymmetry parameters, e.g. sfermion masses and mixing outside of the 
stop sector. Furthermore, no CP-violating phase for the gaugino parameter Mi was included 
in this analysis, which might have interesting effects in the neutralino sector. These issues 
will be studied in future work. 

The present study shows that, even in the challenging case of heavy supersymmetric 
scalars, precise cross-relations between collider physics and cosmological processes can be 
established in order to elucidate some of the main unresolved questions in our understanding 
of the universe. 



Appendix: MSSM case study scenarios 

BGEN: Baryogenesis scenario 

The numerical analysis is based on a typical MSSM parameter point characterized by the 
following weak scale values: 

m | 3 = -gg 2 GeV 2 , Mi = 118.8 GeV, 
m ^ = 4330 GeV, M 2 = 225 GeV, 
A t = -1100 GeV, |/x| = 225 GeV, 

M = O.2, 

m Q,U,D,L,Ri,2 = 10 TeV ' tal1 / 9 = 5 > 

m A o = 800 GeV. 

Due to constraints from large one-loop sfermion-neutralino and sfermion- chargino effects to 
the electric dipole moments of the electron and neutron, the sleptons and squarks of the 
first two generations are chosen to be heavy. The masses of the sbottoms and staus are not 
specified since they have no relevance for the baryogenesis scenario. 

The chosen parameters are compatible with a strongly first order electroweak phase 
transition for electroweak baryogenesis, v(T c )/T c > 1 [5,7], generate a sufficiently large 

baryon asymmetry, rj ~ 0.6 x 10~ 10 , and yield a value for the dark matter relic abun- 
dance 2 of f2cDM^ 2 — 0.109, well within the WMAP bounds. Note that the stop-neutralino 
co-annihilation mechanism is effective for the evolution of the dark matter density in this 
scenario. Furthermore, the stop parameters are chosen such that the mass of the light- 
est Higgs boson is m h o = 117.3 GeV, to satisfy the bound from direct searches at LEP 
my > 114.4 GeV [24]. It was checked that the minimum of the scalar potential is color con- 
serving [7]. At tree- level the following masses are obtained for the relevant supersymmetric 
particles: 

m fl = 122.5 GeV, m^> = 106.6 GeV, m^o = 231.2 GeV, m-± = 162.7 GeV, 
m h = 4333 GeV, m £ » = 170.8 GeV, m^o = 297.7 GeV, m-± = 296.2 GeV, (34) 
cos0 t ~ = 0.010. 

2 The relic dark matter density has been computed with the code used in Rcf. [22]. 
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LCC2: Focus-point scenario 

The point LCC2 is chosen as a point with sizeable gaugino-Higgsino mixing, allowing large 
neutralino annihilation cross-sections into vector bosons. It was studied previously in Refs. [27, 
28]. The scenario is defined by mSUGRA parameters at the unification scale, 

m = 3280 GeV, M 1/2 = 300 GeV, A = 0, sign(^) = +, tan/? = 10. (35) 

Since the evolution of parameters in the focus-point region sensitively depends on the top 
quark mass m t , it is fixed to the value m t = 175 GeV. With Isajet 7.69 [45], the weak 
scale soft breaking parameters are obtained as follows: 



m u 3 


= 1969 GeV, 


Mi 


= 123.1 GeV 




= 2710 GeV, 


M 2 


= 237.6 GeV 




= 3240 GeV, 


l/'l 


= 178.6 GeV 


m u 


= 3268 GeV, 


tan j3 


= 10, 


m R 3 


= 3252 GeV, 




= 3242 GeV, 


Q,U,D,L,Ri, 2 


~ 3300 GeV. 







(36) 



For the relevant supersymmetric particle masses, Isajet 7.69 gives 

m,» = 107.7 GeV, m,o = 190.0 GeV, m-± = 159.4 GeV, m~„ = 850 GeV, 
m.o = 166.3 GeV, m,o = 294.4 GeV, m~± = 286.8 GeV. 

Following the procedure in Ref. [28], the relic density is computed with DarkSUSY 4.1 [44] 
to tt CDM h 2 = 0.109. 
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